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A SPATIALLY EXPLICIT MODEL FOR COMPETITION AMONG 
SPECIALISTS AND GENERALISTS IN A HETEROGENEOUS 

ENVIRONMENT 

By N. Lanchier and C. Neuhauser 

University of Minnesota 

Competition is a major force in structuring ecological communi- 
ties. The strength of competition can be measured using the concept 
of a niche. A niche comprises the set of requirements of an organ- 
ism in terms of habitat, environment and functional role. The more 
niches overlap, the stronger competition is. The niche breadth is a 
measure of specialization: the smaller the niche space of an organism, 
the more specialized the organism is. It follows that, everything else 
being equal, generalists tend to be more competitive than specialists. 
In this paper, we compare the outcome of competition among gen- 
eralists and specialists in a spatial versus a nonspatial habitat in a 
heterogeneous environment. Generalists can utilize the entire habi- 
tat, whereas specialists are restricted to their preferred habitat type. 
We find that although competitiveness decreases with specialization, 
specialists are more competitive in a spatial than in a nonspatial 
habitat as patchiness increases. 

1. Introduction. An ecological community is an assembly of populations 
that inhabit the same spatial location at the same time. Understanding the 
forces that shape such communities has been the goal of community ecology. 

Competition was identified early on as an important driver of community 
assembly and is closely linked to the concept of a niche, which was introduced 
by Elton in 1927 [12]. Elton defined niche as "the status of an animal in its 
community" and used it to explain how multiple species can coexist within 
a community. Hutchinson [15] formalized this concept and defined niche 
mathematically as a subset of an n-dimensional hypervolume. Species with 
overlapping niches compete with each other. One of the main insights has 
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been that in order for species to coexist, they cannot occupy identical niches, 
that is, identical subsets of the niche space. Although delineating the exact 
niche of a species is nearly impossible in nature, the concept has proven quite 
useful, particularly in mathematical models where niches can be identified 
more easily. 

Species that coexist are able to divide up the niche space so that each 
species occupies some subset of the niche space by itself. It thus follows that 
communities composed of species with narrower niches can be more species- 
rich than those composed of species with broader niches (see, e.g., [6]). We 
call the former specialists and the latter generalists. 

Much research has been done to understand the roles of generalists and 
specialists in ecological communities. No clear patterns have emerged. Some 
taxa consist of mostly specialists. For instance, about 90% of herbivorous 
insects are highly specialized, feeding on three or fewer plant families [2]. 
Others are dominated by generalists. For instance, 60% of human pathogens 
are zoonotic, that is, they infect both humans and nonhuman animals, and 
often infect hosts from different orders or classes [22]. 

Evolutionary ecologists have used optimal behaviors of individuals to gain 
insights into the question of relative success of specialists versus generalists, 
in order to assess these strategies in an evolutionary context and ultimately 
to predict community structure. Extending the work of Rosenzweig [21], 
Brown [5] found in a theoretical study on habitat selection that in an envi- 
ronment with two habitat types, three scenarios are possible: (a) two special- 
ists, one for each habitat type, (b) one generalist who occupies both habitat 
types and (c) one generalist and one specialist. Scenario (a) occurs when 
there is no cost to habitat selection, whereas scenario (b) occurs when there 
is a cost. Scenario (c) occurs when there is an asymmetry in the two habitat 
types, either in terms of relative abundances or productivities. 

The theoretical work mentioned above relied on nonspatial models and 
thus did not incorporate a way to account for the often very localized interac- 
tions between competitors. Past research on competition has demonstrated 
that the inclusion of a spatial component with local interactions can alter 
predictions based on nonspatial models [19] or offer new insights into con- 
sequences of local interactions [1]. In this paper we employ the multitype 
contact process [18] in a heterogeneous environment to investigate how lo- 
cal, competitive interactions among specialists and generalists affect their 
abilities to coexist. We define the environment as a patchwork of two host 
types on which specialists and generalists feed. Specialists feed on one of the 
two hosts, whereas the generalists can consume either host. We assume that 
the environment does not change over time. 

We begin by defining the environment. To set the spatial configuration of 
the host population, we let L > 1 be an integer, which we will refer to as the 
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space scale, introduce the ci-dimensional box TLl = [—L,L) d and define 7i\ 
and Ti.2 as a partition of Z d given by 

Hi= |J {2Lz + H L ) and H 2 = \J (2Lz + H L ), 

where H± and H2 denote the sets of z G Z d such that 01 + • • • + z& is even 
and odd, respectively. In our biological context, TCi will represent the set of 
sites that are permanently occupied by a host of type i, i = 1, 2. 

Into this host population we introduce two specialist species and one 
generalist species that differ in their abilities to reproduce. We will refer 
to these species collectively as consumers. To describe their evolution, we 
consider a continuous-time Markov process whose state at time t is a function 
£t : Z d — > {0, 1, 2, 3}. A site x G Z d is said to be occupied by a specialist of 
type 1 (resp. 2) if = 1 (resp. 2), and a generalist if = 3, with the 
condition £(x) = denoting the absence of a consumer. The local dynamics 
at site x are described by the following transition rates: 

— ► 1 at rate a l{f(«)=i;xeWi} 1 — ► at rate 1, 

0<||x-z||<fl 

— ► 2 at rate a !{5(z)=2;xeW 2 } 2^0 at rate 1, 

0<||a;-z||<R 

0^3 at rate /? 1{^( Z )=3} 3^0 at rate I, 

0<\\x-z\\<B. 

where ||x|| = sup i=1 2 ^ l^tl- We refer to R as the dispersal range. To under- 
stand the evolution of our stochastic process when d = 2, let us imagine an 
infinite chessboard where each of its squares has length 2L (and thus 2L x 2L 
sites). The black and white squares represent parts of the habitat that are 
permanently occupied by hosts of type 1 and 2, respectively. Particles of 
type 3, that is, the generalists, do not see the colors of the chessboard and, 
in the absence of particles of type 1 and 2, perform a (translation-invariant) 
contact process. That is, each 3 tries to give birth onto each of its neighbor- 
ing sites at rate (5, the birth occurring if and only if the offspring is sent to a 
site that is not already occupied by a consumer, regardless of the type of the 
host present at that site. The dynamics of the specialists are quite different 
since black (resp. white) squares present unsuitable habitat for particles of 
type 2 (resp. 1). Specifically, if the offspring of a 1 (resp. 2) is sent to a white 
(resp. black) square, then the birth is suppressed. 

The mean-field model. Before we describe the behavior of the spatially 
explicit stochastic model, we will look at the mean-field model [9], that is, we 
will pretend that all sites are independent and that the probability of a site 
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being occupied by a consumer of type i depends only on the host type at that 
site, but is otherwise spatially homogeneous. This then results in a system 
of differential equations for the densities of specialists and generalists. Let 
Vij denote the density of hosts of type i associated with a consumer of type 
j, where i = 1, 2 and j = 1, 2, 3. We let 

ui = ~ - Vu - v 13 and u 2 = \ - v 2 2 ~ v 2 3 

denote the densities of unassociated hosts of type 1 and 2, respectively. The 
constant 1/2 corresponds to the limit of the proportions of the two host types 
in a spatial box whose size tends to infinity. The mean-field limit is obtained 
by letting the dispersal range R tend to infinity. To obtain a meaningful 
limit, we need to rescale the parameters a and (5 by the neighborhood size, 
that is, we set 

a = — and f3 = — , 

where vr denotes the size of the dispersal neighborhood, that is, vr = 
|{z:0 < ||x — z\ < R}\. The following system of differential equations de- 
scribes the mean-field behavior: 

vn + auivn, 

V22 + aU 2 V22, 
VIZ + foil (^13 + ^23), 



dvu 
dt 

dV22 

dt 
dvi3 
dt 




Fig. 1. Solution curves of the mean-field model. Picture on the left: a = 5 and b — 2. 
Picture on the right: a — 3 and 6 = 2. In each of the pictures, we start with vu = V22 and 
«13 = V23 ■ 
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—jj- = -V23 + bu 2 (vi3 + V23)- 

A standard linear stability analysis reveals the equilibrium behavior [20]. 
If b < 1 and a < 2, then the trivial equilibrium [v\\ = V22 = V13 = ^23 = 
0) is the only locally stable equilibrium. If a > 2b and a > 2, then there 
is a locally stable equilibrium point with v\\ = v 2 2 > and V13 = V23 = 
(Figure 1; picture on the left). If a < 2b and b > 1, then there is a locally 
stable equilibrium in which v\\ = V22 = and V13 = V23 > (Figure 1; picture 
on the right). Furthermore, if one of the specialists is initially absent (say, 
V22 = at time 0), then there is an additional equilibrium in which vu, V13 
and U23 are positive, provided a > 2 and < b < § ■ This equilibrium can, 
however, be invaded by specialist 2, which then outcompetes the generalist, 
resulting in the equilibrium in which v\\ = V22 > and V13 = V23 = 0. There 
is no locally stable equilibrium in which both specialists and the generalist 
coexist (i.e., all three have positive densities). A necessary condition for such 
an equilibrium is a = 2b > 2, but the resulting equilibrium is not locally 
stable (it is neutrally stable). 



The spatially explicit model. We now return to the spatially explicit 
model with parameters a and (3. Our results show that the inclusion of 
a spatial structure in the form of local interactions can significantly modify 
the outcome of competition among specialists and generalists. While in the 
mean-field model, both the dispersal range R and the space scale L lose 
their meaning, the study of the spatial model reveals that setting L much 
larger than R helps the specialists to compete. 

As explained above, in the absence of specialists, the process for the gen- 
eralists reduces to a translation-invariant contact process with parameter (5. 
In particular, there exists a critical value A c E (0, 00) such that the following 
holds: If P < A c , the process converges in distribution to the "all 0" configu- 
ration. If [3 > A c , there exists a stationary measure 113 that concentrates on 
configurations with infinitely many 3's ([3] or [17], Theorem 2.25). A sim- 
ple coupling argument shows that specialists die out when a < A c . To show 
this, we consider an initial configuration £0 with £,o(x) ^ 3 for all x £ 
Let £t denote the contact process with parameter a and initial configuration 
Co 0*0 = if £o(x) = 0, and Co(^) = 3 otherwise, that is, Co can be deduced 
from £0 by replacing each specialist by a generalist. Both processes can then 
be constructed graphically from the same collection of independent Poisson 
processes [13] in such a way that Ct has more consumers than £t at any time 
t > (see Section 2 for a detailed description of this graphical representa- 
tion). Since Ct (the process for the generalists) dominates Ct (the process for 
the specialists) and Q converges to the "all 0" configuration for a < A c , it 
follows that Ct also converges to the "all 0" configuration for a < A c . Hence 
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specialists die out when a < A c . It is more difficult to show that specialists 
can survive when a > A c (see Theorem 3). 

To investigate the competition between specialists and generalists, we now 
consider the process starting with infinitely many consumers of each type. 
Relying on a new coupling argument, we may run our stochastic process £t 
and the 3-color multitype contact process rjt (in which particles of type 1 
and 2 give birth at rate a and particles of type 3 give birth at rate (3) on 
the same probability space in such a way that, starting from the same initial 
configuration, £j has more 3's and fewer l's and 2's than rjt. According to 
Theorem 1 in [18] the process rjt converges in distribution to 113, provided 
(3 > a > A c . The above-mentioned coupling then implies that the generalists 
outcompete the specialists. An analysis of the dual process will show that 
the result still holds when a = (3 > A c . That is, while coexistence occurs for 
the multitype contact process in a homogeneous environment when a = j3 
and d > 3 (see [18]), the process in this heterogeneous environment gives the 
generalists an advantage over specialists in any dimension. This is summa- 
rized in the following theorem where denotes weak convergence, and 
fj,3 is the upper invariant measure of the basic contact process introduced 
above. 

Theorem 1. Assume that L and R are fixed and that (3>a and (3 > A c . 
If £o is translation-invariant and P(£$(x) = 3) > 0, then £t => fJ>3- 

This result corresponds to the case b > a and b > 1 in the mean-field model 
and tells us that the fragmentation of the environment allows the generalists 
to outcompete the specialists in both the spatial and the nonspatial models, 
provided the birth rate of the generalists is at least as high as that of the 
specialists. 

The next step is to investigate the behavior of the spatially explicit model 
when the specialists have a higher birth rate than the generalists. In the 
mean-field model, we have seen that the condition a > 2b is required for the 
specialists to invade the generalists in their equilibrium. The behavior of the 
stochastic process, however, is different. Let us choose a space scale L much 
larger than the range of interactions R so that well inside black squares, 
specialists of type 1 and generalists behave nearly like a multitype contact 
process. In particular, specialists locally outcompete generalists, provided 
a > (3. The survival of specialists, however, is not clear since the contact 
process on a finite set (here the 2L x 2L square when d = 2) always converges 
to its absorbing state, the "all 0" configuration. To show that the specialists 
actually survive on the infinite chessboard (with large squares), we will prove 
that, with probability close to 1, the l's can invade the four diagonally 
adjacent black squares before going extinct. The trickiest point of the proof 
will be to show that, with high probability, a specialist of type 1 in the middle 
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of a black square will produce four "invasion paths" to bring its offspring 
within a finite distance (independent of L) of each of the four corners of 
the square. Since white squares are forbidden for the l's, the invasion paths 
will have to be contained in the black square. Comparing the set of black 
squares occupied by (a significant quantity of) l's with the set of wet sites in 
a one-dependent oriented percolation process will then yield the following: 

Theorem 2. Assume that R is fixed and that initially there are infinitely 
many l's and 3's, but no 2's. If a > /3 > A c and d>2, then l's and 3's 
coexist, provided L is sufficiently large. 

Since specialists of type 1 and 2 do not interfere with each other, the proof 
of Theorem 2 also implies the following: 

Theorem 3. Assume that R is fixed and d>2. If a > A c and a > j3, 
and initially there are infinitely many l's and 2's, then l's and 2's coexist, 
provided the space scale L is sufficiently large. 

Unfortunately, Theorem 3 does not tell us if all three consumers coexist 
when d>2. Numerical simulations on a 200 x 200 torus indicate that the 
generalists persist for a very long time along the boundaries of the squares of 
the chessboard, namely where the density of specialists is low, but proving or 
disproving that they can percolate on the infinite lattice would be difficult. 

Comparison of the spatially explicit and the mean-field models. In the 
mean-field model, that is, in the absence of spatial structure, the condition 
for the specialists to invade the generalists is given by a > 2b. The factor 2 
comes from the fact that only half of the habitat is available for each spe- 
cialist. In the spatially explicit model, the dimension d as well as the space 
scale L and the dispersal range R can profoundly modify the outcome of the 
competition. In the one-dimensional case, the behavior of the spatial model 
is quite easy to predict when R < 2L. Since all the offspring of a specialist 
are confined to their parents' segment of length 2L < oo and are unable to 
reach other segments (because R < 2L), the specialists eventually go extinct 
so that the only nontrivial stationary distribution of the process is ^3. To 
describe the behavior of the process when d>2, let a c denote the infimum 
of a's such that the specialists can survive on the chessboard. Theorems 
1 and 3 tell us that if the birth rate of the generalists is equal to (3, then 
(3 < a c < co. Moreover, a c — > (3 as the space scale L tends to infinity. That 
is, taking large L makes the specialists much more competitive in the spa- 
tially explicit model than in the mean-field model. The right-hand picture 
in Figure 2 illustrates the effects of L on the outcome of the competition. 
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Fig. 2. Simulation results for the spatially explicit model with parameter p = 2 on f/ie 
200 x 200 square with periodic boundary conditions. In both pictures, the vertical axis 
denotes the common birth rate of specialists. The curves represent the smallest a such 
that the density of specialists is greater than the density of generalists after 100 units 
of time. The simulation results suggest that when the range of the interactions of the 
symbionts is larger than the length of the squares, the competitiveness of the specialists 
reduces significantly, just as in the mean-field model. The competitiveness of the specialists 
increases when the range interactions is short and the habitat is coarse-grained as shown 
in Theorem 2. 

Another interesting corollary of Theorem 3 is that, in the absence of gen- 
eralists, A c < a c < oo and a c — > A c as L tends to infinity. In conclusion, our 
(multitype) contact process on the chessboard behaves nearly like the corre- 
sponding (multitype) contact process in a homogeneous environment when 
the size of the squares is large, while their respective nonspatial versions 
differ significantly. 

2. Duality. Proof of Theorem 1. We saw in the Introduction that the 
behavior of the process when a < (3 is a straightforward consequence of 
Theorem 2 of [18], so we will now focus on the case a > (3. The key to the 
proof of Theorem 1 is duality. To define the dual process of our heterogeneous 
multitype contact process, we start by introducing a graphical representation 
from which the process with parameters a and (3, where a > j3, can be 
constructed. 

The graphical representation when a> j3. For any x, z G Z d , < \\x — 
z\\ < R, let {T£ ,z : n > 1} and {U% : n > 1} denote the arrival times of inde- 
pendent Poisson processes with parameters a and 1, respectively. At times 
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T£' z , we draw an arrow from site x to site z, toss a coin with success prob- 
ability (a — (5) I a and, if there is a success, label the arrow with an s (for 
specialist). Moreover, if x £ Hi and z £ TCj with i^j, then we label the 
arrow with a g (for generalist). This can be interpreted as follows. If, at 
time T£' z , the site x is occupied by a consumer and z is empty, then the 
consumer at site x, regardless of type, gives birth to a consumer of the same 
type through this arrow, unless the arrow is labeled with an s or a g. Only 
specialists (resp. generalists) can give birth if the label is s (resp. g). We 
stipulate that arrows labeled with both s and g cannot be used by either 
consumer. To complete the construction, we put a "x" at (x, U%) to indicate 
that a consumer, regardless of its type, dies. 

The dual process. The dual process in the heterogeneous environment 
is defined as for the process in a homogeneous environment. We say that 
there is a path from (x,s) to (z, T), where < s < T, if there is a sequence 
of times sq = s < s\ < ■ ■ ■ < s n +i = T and a sequence of spatial locations 
xq = x,x±, . . . ,x n = z such that the following two conditions hold: 

1. for i = 1, 2, ... ,n, there is an arrow from a?j_i to Xi at time Sj, and 

2. for i = 0, 1, . . . , ra, the vertical segments {x^ x (sj,Sj+i) do not contain 
any x's. 

A dual path from (x, T) to (z, T — s) then refers to a path from (z, T — s) to 
(x,T). The dual process starting at (x,T) is defined by setting 

|(*,T) = { ze z d : there is a dual path from (x, T) to {z, T - *)}. 

The dual process is naturally defined only for < s < T but it is convenient 
to assume that the Poisson processes involved in the graphical representation 
are defined for negative times so that the dual process can be defined for all 
s>0. 

We introduce the dual process in order to deduce the state of site x at 
time T from the configuration at earlier times. The aim is to understand 
the ancestry of (x,T) by working backwards in time. As for the multitype 

contact process, at any time s > 0, the elements of £g X ' , called ancestors, 
can be arranged according to the order in which they determine the type of 
(x,T). This hierarchy can be deduced from the tree structure of 

r = {(H x > t) ,s):0<s<T}, 

and we refer to [18], Section 2, and [11], Section 3, for a detailed description. 
To have a rigorous definition of the hierarchy of ancestors, we introduce a 
function cj) s that, at any time s > 0, maps the dual process into the set of 
sequences with integer values S equipped with the usual lexicographic order 
<S- We recall that for any v,w G 5, 

v <C w if and only if Vf. = Wk for k = 1, 2, . . . , ri — 1 and v n < w n 
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for some integer n > 1. To facilitate the writing of 4> s , we identify any v £ S 
such that 

v n 7^ and f fc = for > n + 1 

with the vector (v±, V2, ■ ■ ■ ,v n ). The function <f> s is defined inductively in the 
following way. First, we go down the graphical representation starting at 
(x,T), denote by T — s x the first time a death mark is encountered at site 
x and let 4> s (x) = for any < s < s x . Assume that a site z is added to the 
dual process at time T — u and denote by T — s z the first time we encounter a 
death mark at site z by going backward in time starting from (z, T — u). We 
then go back up from (z, T — s z ) to (z, T — u) and, each time we encounter 
the tip of an arrow, denote its space-time location by {z k ,T — Sk), with 
s z > si > • • • > s n > u. For k = 1, 2, . . . , n, denote by T — s Zk the first time a 
death mark is encountered at site z^. We then set 

<Ps(zk) = (u,k) where u = 4> Sk (z) fois k <s<s Zk . 

See Figure 3, left-hand picture. This naturally introduces a hierarchy in the 

"(x T) 

set of ancestors where, given two sites y and z belonging to Q ' , site y 
comes before site z in the hierarchy if and only if (f) s (y) <C (f> s (z). The spatial 

location of the first ancestor at time t — s is denoted by ^'^(1). Later on, 
this ancestor will be called the distinguished particle. 

The space-time point (x, T) is said to live forever if ^ :C ' T ^ 7^ for all s > 0. 
Provided (x,T) lives forever, the distinguished particle has a nice feature, 
namely that its path can be broken up at renewal points into independent 
and identically distributed pieces. Renewal points are defined as follows. If 
at time s the distinguished particle jumps to a site y such that (y, s) lives 
forever, then (y,s) is a renewal point. Let (Sk,Tk) denote the kth renewal 
point, Xi the spatial displacement between consecutive renewal points and 
Tj the corresponding temporal displacement, that is, 

k k 

Sk=X + and T k = Ti - 

i=l i=X 

See Figure 3, right-hand picture. 

Proposition 2.1 ([18]). Conditioned on the event that (x,T) lives for- 
ever, the random vectors {(-Xj, Tj)}j>i are independent and identically dis- 
tributed. Moreover, the tail distributions of Xi andri have exponential bounds, 
that is, 

P{\\Xi\\ >t) <Ce~^ and P(n > t) < Ce^ 
for appropriate C < 00 and 7 > 0. 
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u. i) 



(2. i) 



(2.1.21 



(1.2. I) 



(I. 1,1) 



(2. I. IJ 



(1.1,1) 



I - I - J I 



(J. 1. 2. t) 



Fig. 3. Tree structure of the dual process. Left-hand picture: Construction of the function 
4> s and description of the ancestor hierarchy at time 0. Right-hand picture: Spatial and 
temporal locations of the renewal points. In both pictures, the bold lines refer to the path 
of the distinguished particle. 



Before proving Theorem 1, we need to describe an algorithm that allows 
us to deduce the type of (x,T) from the (topological) structure of T and the 
hierarchy of the ancestors. We look for the first ancestor in the hierarchy that 
does not land on an empty site at time 0. If this ancestor lands on a specialist 
(resp. generalist) at time 0, then (x,T) will be occupied by a specialist (resp. 
generalist) of the same type unless the path crosses a g-arrow (resp. s-arrow) 
on its way up to (x,T). If it lands on a specialist (resp. generalist) and the 
path crosses a g-arrow (resp. s-arrow), this specialist (resp. generalist) can 
block some other ancestors from determining the type of (x,T) so we will 
need to remove these ancestors from the hierarchy. To do this, we follow the 
path on its way up to (x, T) until we first encounter a g-arrow (resp. s-arrow), 
then remove from the hierarchy all the ancestors of the dual process starting 
at the point where this arrow is attached, repeat the same algorithm with the 
first remaining ancestor of the hierarchy at time 0, and so on. The algorithm 
halts whenever we find an ancestor that successfully reaches (x,T). If none 
of the ancestors succeed, (x,T) will be empty. 



Proof of Theorem 1. It remains to prove Theorem 1 for a = (3. The 
proof is similar to the proof of Theorem 1 in [18]. We begin by assuming 
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that the dual process starting at (x,T) lives forever, an event with positive 
probability since a > A c . We will prove that P(£t(x) = 3) — ► 1 as T — > oo on 
the set where (x,T) lives forever. Without loss of generality, we can suppose 
that x £ H\. Proposition 2.1, together with the irreducibility of the random 
walk associated with the renewal points, implies that there is a subscript 
k almost surely finite such that € Ti.2- In particular, since what happens 
before and after each renewal is independent, we can inductively extract a 
subsequence of renewal points as follows. We set k$ = and, for any i>l, 



By continuity of the path, the distinguished particle will cross, on its way 
up to (x,T), at least one arrow that crosses from TC\ to H2 (or vice versa) 
between time T — T^. and time T — T) c ._ 1 . It follows from our graphical 
representation that such an arrow is a g-arrow which will block the specialists 
from determining the type of (x,T). 

To prove that P(£t(x) = 3) — > 1 as T — > 00 [provided (x,T) lives forever], 
we will construct a sequence of ancestors that are good candidates to bring 
a generalist to (x,T). We start the dual process at (x,T) and determine the 
ordered set of ancestors after T units of time by going backward in time. The 
first member of the sequence, denoted by ??r(l), is the distinguished particle. 
We follow the path this ancestor takes to determine the type of (x,T) by 
going forward in time until the first time we cross a g-arrow. We then discard 
all the ancestors that cross this g-arrow on their way up to (x,T). The first 
ancestor that is left after discarding these ancestors, denoted by ?7t(2), is the 
second member of our sequence. Since, before discarding, the tree starting 
at each renewal point is linearly growing in time (in the sense described by 
the Shape Theorem in [17], page 128), we can find, by choosing t sufficiently 
large, an integer i > 1 such that, even after discarding, the dual starting at 
(SktjTfy) lives forever. This tells us that the path the second member will 
take to determine the type of (x,T) goes through the kith, renewal point. 
The definition of the sequence of times , i > 1 , implies that the second 
member crosses at least ki g-arrows on its way up to (x,T). As previously, 
we discard all the ancestors that cross the first g-arrow encountered by 
the second member on its way up to (x,T), define the third member of 
the sequence as the first remaining ancestor, and so on, until we run out of 
ancestors. Since different ancestors can occupy the same site, we now extract 
a further subsequence so that, at time 0, all the candidates are different. We 
start with 77^(1) and discard all the members of {r]T(k) :k>2}, that occupy 
the same site as r/r(l)- We then take the next ancestor that is left, and so 
on. 



let 




min{A; > : Sk G H2] 
min{A; > : Sk £ Hi} 



if i is odd, 
if i is even. 
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We now denote by rjr the set of members of the subsequence we have 
just defined. With the possible exception which can occur if it is the last 
member of the subsequence [since it may happen that this member does not 
cross any g-arrow on its way up to (x,T)], if a site z £ t]t is occupied by 
a specialist at time 0, then this specialist will be blocked on its way up to 
(x,T) by a g-arrow. On the other hand, if z is occupied by a generalist, this 
generalist will be able to bring its offspring up to (x,T) since there are no 
s-arrows in our graphical representation (recall a = (3). In particular, if we 
denote by 63 the set of sites occupied at time by a particle of type 3, we 
only need to prove that 

P(VT n G 3 / 0) -» 1 asT^oo. 

Since the tree growing out of (x, T) is linearly growing in time, for any e > 
and M > 0, there exists a time to > such that 

P(card??T < M) < e VT > t . 

To conclude, we apply Lemma 9.14 from [14] which tells us that if £0 is 
translation- invariant and P(£o{x) = 3) > then, given e > 0, there is an 
M > such that if cardr/y > M, then 

P(r) T n e 3 = 0) < e. 

This, together with the previous inequality, completes the proof of Theorem 
1. 

3. Proof of Theorem 2. This section is devoted to the proof of Theorem 
2, which describes the behavior of the stochastic process when a > (3 > A c 
and the space scale L is large. To avoid cumbersome notation, we will prove 
the result only when d = 2, but would like to point out that when d > 3, 
each step of the proof is the same. The proof relies on a rescaling argument 
(see, e.g., [4, 8]). The basic idea is to show that, for any e > 0, members of 
the family of processes under consideration, when viewed on suitable length 
and time scales, dominate a one-dependent, oriented site percolation process 
in which sites are open with probability at least 1 — e. 

Let B = TCl = [—L,L) 2 . To prepare for a multiscale argument, we define 

a sequence of spatial boxes J n = [—L/n,L/n) 2 , n = 1,2, For any z = 

{z\,Z'i) in Z 2 and any positive integer n G N*, we set 

&(z) = (2Lz 1 ,2Lz 2 ), B(z) = $(z) + B and J n (z) = &(z) + J n . 

Moreover, we tile J&{z) with L 0,1 x L 0,1 squares by setting 

7t(w) = (L°- 1 w 1 ,L 0A w 2 ), D = [-L 0A /2,L 0A /2) 2 , 

D{w) = tt(w) + D, I z ={weZ 2 : D{w) C J±(z)}. 

To prove that specialists and generalists coexist, we start by introducing the 
following two collections of good events. We set T = L 1 ' 5 . 
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1. For z\ and z 2 both even or both odd, we will say that J±(z) is s-good if 
Ji(z) is void of 3's and has at least one particle of type 1 in each of the 
squares D(w), w € I z . For z\ and z 2 both even for even k and z\ and Z2 
both odd for odd k, we will say that (z, k) is s-occupied if the spatial box 
Ji(z) is s-good at time kT. 

2. For z\ even and z 2 odd or z\ odd and z 2 even, we will say that J±{z) is 
g-good if J4(z) is void of l's and has at least one 3 in each of the squares 
D(w), w € I z - For z\ even and z 2 odd for even k and zi odd and z 2 even 
for odd k, we will say that (z,k) is g-occupied if the spatial box J4(z) is 
g-good at time kT. 

The main ingredient needed to prove that the specialists survive is the fol- 
lowing: 

Proposition 3.1. Assume that a > j3 > X c and J4 is s-good at time 0. 
For any e > 0, the space scale parameter L can then be chosen in such a way 
that 

P(J&(1, 1) is s-good at time T) > 1 — e. 

The idea is that if some square B(z) C TLi is occupied by a cluster of 
specialists at a given time, then this cluster will invade the four diagonal 
squares with probability close to 1, and this will occur in less than T units 
of time. It then follows from well-known results about oriented site percola- 
tion (see, e.g., [7, 10]) that the s-occupied sites percolate provided e > is 
sufficiently small, implying that the specialists survive. To prove that both 
consumers coexist, we still need to show that the g-occupied sites percolate 
as well. More precisely, we have, analogously: 

Proposition 3.2. Assume that (3 > X c and J<i(Q, 1) is g-good at time 
0. For any e > 0, the space scale parameter L can then be chosen in such a 
way that 

P(J4(1,0) is g-good at time T) > 1 — e. 

We illustrate the strategy of Proposition 3.1 in Figure 4. Let z = (1, 1) and 
x G Ja(z). We will then construct a dual path At starting at Aq = (x, T) that 
lands at time on the target set J4 in the heart of the cluster of specialists. 
This selected path will have to block 3's, but not l's, from determining the 
type of (x,T). In particular, the important step will be to move At from 
Ja(z) to J4 without escaping from the set Hi, by crossing the corner at the 
point Lz = (L,L). The proof of Proposition 3.2 is similar. 
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Construction of the selected path. The construction of the selected path 
A t relies on the repositioning algorithm described in [11], supplemented with 
a restart argument. To define At, we will first construct a sequence of auxil- 
iary processes for k = 0, 1, ... , denoted by {£,t}t>o, called the kth test path, 
with ^ taking values in 1? , and then show in a series of lemmas that we 
can use the test paths to construct At. We will say that a renewal point is 
associated with an s-arrow if the first arrow a particle crosses, starting at 
this renewal point and moving up the graphical representation, is an s-arrow 
(see Figure 5, left-hand picture). Furthermore, let A denote the straight line 
going through the origin (0, 0) and the point Lz, and denote by dist(-, ■) the 
Euclidean distance. We also let 5 > be a constant, to be fixed later. 

The test paths are defined inductively (see Figure 5, right-hand picture). 
To construct the first test path {Ct}t>o, we define a sequence of stopping 
times aij, j = 0, 1, 2, . . . , with £71,0 = 0. The process £,} starts at £q = x an d 
follows the path of the first ancestor starting at (x,T) until the first time, 
denoted by o"i,i, that it jumps to a renewal point associated with an s-arrow. 
At time o~i t i, we either leave f^ 1 where it is at that time or we reposition 
To determine whether and where to reposition we denote the location 
of the second ancestor in the hierarchy at time o~i\, if it exists, by B ai 1 . If 
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Fig. 5. Left-hand picture: Picture of the renewal points. The first renewal point (Si,Ti) 
is not associated with an s-arrow. The second one, namely (E>2,T2), is associated with an 
s-arrow. In particular, ai,i = T%. Right-hand picture: Picture of the first test path (solid 
path). The bold straight lines refer to the jump process St ■ The points represent the renewal 
points associated with an s-arrow. The test path £| is repositioned at closed points and left 
where it was at that time at open points. 



B ai 1 exists and the dual process starting at (B ai 1 ,T — <7i,i) lives forever, 
then we have the following three alternatives: 

1. Ifdist(|^ 1 ,A)>5anddist( J B (T11 ,A)<dist(|^ 11 ,A),thenweset|^ il+ = 

2. If dist(|£ M , A) < 5 and dist(5 t7l A , 0) < dist(|^ x , 0), then we set ^ 1+ = 

3. Otherwise, we set ^ 1+ = ^ 1 . 

If B ai 1 does not exist, or the dual process starting at (B ai l ,T — ci,i) does 
not live forever, then we set taii+ = Wi i • 

To define the next stopping time (J\ t 2, we start a new dual process at 
(£^1 i + ' T ~ an d follow the path of the first ancestor until the first time 
it jumps to a renewal point that is associated with an s-arrow. Denote this 
time by (J\,2- [Note that we followed the path of the first ancestor of the dual 
process starting at 1+ ,T — cri t \) for cr lj2 — a^i units of time, as shown 
in Figure 5, right-hand picture.] We apply the repositioning algorithm again 
as described above. 
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Fig. 6. Picture of the first three test paths. The selected path A t , in bold lines, is equal to 
the third test path, which is the first successful path. As in Figure 5, right-hand picture, the 
points represent the renewal points associated with an s-arrow. The repositioning algorithm 
is applied successfully at filled points. 

Into the first test path, we embed a jump process Sj, t> 0, as follows. 
The process Sj stays put except at the stopping times cry, j = 0,1,..., that 
is, Sf = ^ . + for cry <t< 0"y +1 . The effect of the repositioning algorithm 

is that, whenever possible, we move the jump process Sj closer to while 

staying within a finite distance 6 of A. 

We set Tq = and let T/ denote the first time jumps inside Jg. We call 

£t a successful path if a particle, starting at {^Lj,T — T{) and going up to 

1 1 

(x,T) by taking the path we have just defined, does not cross any g-arrow. 
If the 1st path is not successful, we define a 2nd path £f as follows. For the 
first T-f units of time, the test path follows the path of the distinguished 
particle starting at £q = x, and from time T/ on moves according to the 
repositioning algorithm. Into this process, we embed the jump process Sf 
and denote by T 2 J the first time Sf jumps inside J%. We repeat the procedure 
until the first successful path is found, and call this path A t (see Figure 6). 

Existence of a successful path. Since the kth test path $ follows the 
distinguished particle starting at [x, T) from time to time Tj^_ 1 , the first 

step is to prove that, with probability close to 1, the first ancestor ^ x,T \l) 
does not escape from the box Ji(z) by time T = L 1,5 . This allows us to 
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locate the path when we start applying the repositioning algorithm. It also 
implies that our fcth path does not cross any g-arrows between time and 
time Tj?_ 1 and is therefore a successful path. 

Lemma 3.3. Assume that a> j3 > X c . If x G J/l(z) and T = L 1 * 5 , then 

P(£T' T) (1) i J 2 (z) for t G [0,T]) < C t exp(- 7l L - 2 ) 
for appropriate C\ < oo and 71 > 0. 

Proof. Since our graphical representation is invariant under translation 
of the vector 2Lz, we can assume that x G J4 and work with J 2 instead of 
Ji(z). We will prove that, with probability close to 1, (i) the renewal points 
do not escape from the box x + Js and (ii) between two renewal points, the 
spatial displacements of the distinguished particle do not exceed L/8, an 
event that we define as 

G k = {\S k - e?' T) (l)| < L/8 for T k < t < T k+1 }. 

Let mi = Eti and denote by St the location of the renewal point S n at time 
t, and by Nt the number of renewal points between time and time T. We 
then have 

P(£l X ' T \l) i J2 for some t < T) 

<P(N T >2T/mi) 

+ P(G% for some < k < N T ; N T < 2T/m 1 ) 

+ P(S t -x(/Jg for some < t < T; N T < 2T/m x ). 

By using large deviation estimates, we get 

P{N T > 2T/mi) < C 2 exp(- 72 T) 

for appropriate C 2 < 00 and 72 > 0. Furthermore, the proof of Proposition 
2.1 implies that 

P{G% for some < k < N T ; N T < 2T/m 1 ) 

< P(G% for some 1 < k < 2T/m x ) 

< 2T/m 1 P{G\) < 2T/?niCexp(-7L/8). 

Finally, by letting denote the 1st coordinate of St and using the reflec- 
tion principle, we get 

P(S t -x(/J 8 for some < t < T; N T < 2T/m 1 ) 

< 2P(|5 t (1) - 5^1 > L/8 for some < t < T; N T < 2T/m l ) 
<4P(|5« -5«|>L/8). 
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Now, Chebyshev's inequality implies that, for any 9 > 0, 

2T/m 1 

P(\sS /mi -S^\>L/8)<e- &L / s J] Eexp(^f)) 

i=l 

/ 2L L5 , \ 

<exp -— + log 0(0 , 

V 8 mi / 

where 0(0) is the moment-generating function of X^ . Since EX[ 1] = and 

oo, we have log 4>(9) < C36 2 for some C3 < 00 and for sufficiently 
small 9. By taking 9 = L~ ' 8 in the last expression, we can conclude that 

P(S t -x(£J 8 for some < t < T; N T < 2T/m 1 ) < 4exp(-L°' 2 /10) 

for sufficiently large L. Putting things together, we finally get 

Piii^i 1 ) i J 2 for some t < T) 

< C 2 exp(-7 2 L L5 ) + 2L L5 /™iCexp(- 7 L/8) + 4exp(-L a2 /10). 
This completes the proof of the lemma. □ 

The second step in proving the existence of a successful path is to inves- 
tigate the behavior of the process ^J 71 between time Z^_ 1 and time T^, that 
is, when the repositioning algorithm is applied. The aim is to prove that 
there is a positive probability that, between these times, the mth test path 
does not cross any g-arrows on its way to J$. Since this will only happen 
if the path crosses the corner as shown in Figure 4, we will call this event 
to successfully cross the corner. Let t be a positive integer to be fixed later 
and, for any k > 0, set 

a k = Lz + 2 k (£,0), b k = Lz + 2 k (0,l), 

c k = Lz + 2 k (£,2£), d k = Lz + 2 k (2ej). 

Moreover, we consider the sequence of segments 

I k = [3a k /8 + 5b k /8, 5a k /8 + 3b k /8], 

as well as the rapidly increasing sequence of squares given by D k = (a k ,b k , c k ,d k ). 
That is, I k is the segment with endpoints 30^/8 + 56^/8 and 50^/8 + 36^/8, 
and D k the square with corners a k , b k , c k and d k . We can assume that £ and 
L are chosen so that there exists K* with L = 2 K * £. Ik* is then contained 
in the diagonal line connecting the points Lz + (2L, 0) and Lz + (0, 2L) (see 
Figure 7 for a picture). Assume that the mth test path starts at x £ Ik+i- 
We will prove that, with a positive probability (uniform in k), we can bring 
our test path from x to Iq without escaping from the decreasing sequence 
of squares D k , D k _i, . . . , Dq. We say that the repositioning algorithm was 
applied successfully if the new path is chosen. Let Mt denote the number of 
times the repositioning algorithm has been applied successfully by time t. 




FlG. 7. The sequence of squares. 



Lemma 3.4. Assume that a > J3 > A c . There exists 74 > so that 
P{M t < l4 t)<C 5 exp(-j 5 t) 
for appropriate C5 < 00 and 75 > 0. 

Proof. This is the same as the proof of Lemma 3.4 in [11]. □ 

In Lemma 3.5, we will show that, with high probability, a test path can 
move from I\ to Iq without leaving Dq. A corollary to this lemma is that, 
with probability close to 1, a test path can move from to Ik-i without 
leaving Dk-i- Applying this lemma repeatedly will then show that, with 
positive probability, a test path can successfully cross the corner starting 
anywhere in J2{z). 

Lemma 3.5. If a > (3> A c and the mth test path starts within distance 
of I\ , then 



P(the test path escapes from Do before getting within £/16 of Iq) 
<C 6 exp(- 7f /- 4 ) 
for sufficiently large I and appropriate C§<eo and 76 > 0. 
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Proof. To lighten the notation, we will write at instead of a m ^. Let 
9\ (resp. 82) denote the first time dist(S'f l ,A) < 5 [resp. crosses either 
segment (aodo) or (60C0)], and 7ta denote the orthogonal projection on the 
straight line A. It follows from the repositioning algorithm that there exists 
777-2 > such that, on the set where dist(S™ , A) > S, 

E[dist(S™, A) -dist(S™_ i; A)|^ fe _J < -m 2 

where T ( , h _ 1 = &(S™ , . . . , S™ ). This, together with Lemma 3.4, implies 
that we can find constants C 7 , Cg, 78 G (0,oo) such that 

P(01 > CV or 0i > 2 ) < C 8 exp(- 78 £). 

Furthermore, using the notation introduced in Lemma 3.3, we get 

P(dist(4 m , (ai&i)) > £/2 for some t < 0i) 

< P(0i > C 7 £ or 0i > 6 2 ) + P{N ei > 2C 7 £/mi;9 1 < C 7 £) 

+ P(the largest spatial displacement is greater than £/16; 

N 01 < 2C 7 £/m x ) 

+ P(dist(Sf, (0161)) > e/4 for some t < 9 1 ;N Bl < 2C 7 £/m 1 ). 

Since ir^,(S^) — 7r A(5 , ™_ 1 ) has mean (and finite variance) before time 0i, 
the arguments of the proof of Lemma 3.3 (with = £ -0,6 ) together with the 
exponential bound on 0i imply that 

P(dist(C, (ai&i)) > 1/2 for some t < 0i) 

< C$ exp(-7 8 £) + C 2 exp(-j 2 C 7 £) 

+ 2C 7 £/m l C exp(- 7 £/16) + 2 exp(-£ a4 /8). 

Now that £|™ is close to A, the next step is to bring it to the target region 
Iq without escaping from the square Dq. Let 03 > 0i (resp. 04 > 0i) denote 
the first time S™ crosses the straight line (ao^o) [resp. (codo)], and let K 
denote the number of times the repositioning algorithm has been applied 
between time 0i and time 03. First of all, since S™ has a drift toward A 
when dist(S™ , A) > 5, there is a positive constant p such that 

P(cardK e [0i , 3 ] : dist(£™ , A) < 6) < pK) < C 9 exp(- l9 £) 

for appropriate C9 < 00 and 79 > 0. That is, with probability close to 1, the 
fraction of renewal points associated with s-arrows that are within distance 
5 from A is at least p. Furthermore, there exists 7773 > such that, on the 
set where dist(5'^_ i , A) < 5, 

E[dist(S™ 0) - dist(S™, 0)1^] < -m 3 
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which, together with the previous estimate, implies that there is a C\o < oo 
such that 

P(9 3 - 9 1 > C 10 £ or 9 3 > 4 ) < C u exp(- 7 n4) 

for appropriate Cu < oo and 711 > 0. The same decomposition as above 
then implies that 

P(dist(C, A) > £/8 for some 6 1 < t < 63) < C 12 exp(--/ 12 £) 

for appropriate C\ 2 < 00 and 712 > 0. In conclusion, there exist Cq, 76 G 
(0, 00) such that 

P(the test path escapes from Do before getting within £/16 of Io) 

< P(dist(4 m , (0161)) > £/2 for some t<0 1 ) + P{9 1 > C 7 £ or 1 > 6 2 ) 
+ P(dist(|[™, A) > 1/8 for some 9 1 < t < 63) 
+ P(8 3 - 0i > C w £ or 6 3 > 4 ) 
<C 6 exp(- 7f /- 4 ). 
This completes the proof of the lemma. □ 

Relying on Lemma 3.5, we can now prove that there is a strictly positive 
probability that, starting at x G J 2 (z), the mth test path reaches J$ without 
crossing any g-arrow. This, together with independence of temporally sepa- 
rated paths, will imply that a geometric number of trials suffices to find a 
successful path. 

Lemma 3.6. Assume that a > (3 > X c . There exists pi> such that 
Pithe mth test path successfully crosses the corner) >p\ 
for all sufficiently large L. 

Proof. This is a three-step process in which we will prove that, with 
positive probability, we can bring our path from J 2 {z) to Iq, then from Iq 
to Io — (1,1), and finally from Iq — (£,£) to J% without crossing any g-arrows 
(see Figures 4 and 7). To estimate our first good event, we apply Lemma 3.5 
repeatedly to obtain 

P( starting within L/16 of Ik*, 

the path crosses a g-arrow before getting within £/16 of Iq) 

00 

< C 6 J2 exp(-7 6 (2 fc £) a4 ) < 2C 6 exp(- 7f /- 4 ) 

fc=0 
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for sufficiently large t. Note that the previous estimate holds regardless of the 
value of L. Moreover, it tells us that we can fix (. so that our first good event 
occurs with probability p2 > 0. The same reasoning that led to the estimate 
of the first good event can be used to estimate the third good event, since 
it is symmetric to our first good event. We reverse time in Lemma 3.5 and 
find that if a test path starts within £/16 of Iq, then the probability that 
the test path escapes from Dq before getting within £/8 of I\ is less than 
Cq exp(— 76^ ' 4 )- Iterating as before, it follows that, with high probability, 
we can move the particle inside Jg without crossing any g-arrows. Now that 
I is fixed, there is a positive probability ps(l) > that the mth test path 
goes from Iq to Jo — (£, £) without crossing any g-arrow, since the interaction 
neighborhood we chose has (at least) the eight nearest neighbors. Finally, to 
be able to start from anywhere in J2(z), we use the repositioning algorithm 
to move the path first from x S Ji(z) toward A and then toward Ik*- Using 
similar estimates as in Lemma 3.5, it follows that 

P(the test path escapes from J\{z) before getting within L/16 of Ik*) 
<C 6 exp(- 76 L ' 4 ) 

for sufficiently large L and appropriate constants Cq, 76 S (0, oo). In conclu- 
sion, the lemma follows by taking p\ = P2P3(£) > 0. □ 

With Lemmas 3.3 and 3.6 in hand, we are now ready to prove that, with 
high probability, the first successful path reaches Jg by time T — vL. More 
precisely, we have the following: 

Lemma 3.7. Let T = L 1,5 and N = mm{k > is a successful path}. 
It then follows that 

P(T J N >T-^L)< C 13 exp(-7 13 L - 2 ) 

for a> (3 > A c and appropriate C13 < 00 and 713 > 0. 

Proof. The first step is to prove that the temporal displacement Tj! +1 — 
is small enough such that, by taking sufficiently large L, the number of 
test paths run by time T can be made arbitrarily large. The arguments used 
in the proof of Lemma 3.5 to estimate the stopping times 9\ and #3 imply 
that there exists C14 < 00 such that 

P{Ti +1 - T£ > C U L) < C 15 exp(-7 15 L) 

for appropriate C15 < 00 and 715 > 0. In particular, since the increments 
T^ +l — are independent, it remains to be shown that a geometrical num- 
ber of trials suffices to find a successful path. Let Vl denote the event that 
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the distinguished particle starting at (x,T) does not cross any g-arrow by 
time T. From Lemma 3.3 we know that 

P(Q) > l-Ciexp(-7iL a2 ) 

and from Lemma 3.6 that 

P(the kth test path does not 

cross any g-arrow by time T^|0;T^_ 1 < T) > p±, 

regardless of the values of k and L. This implies that 

P(T£ > T - v 7 !) < P{T^ >T-^L;N> L 0A ) 

+ P(T N >T-VZ;N< L oa ; Q) 

+ P{T N >T-VZ;N< L oa ; Q c ) 

< (1 - pi) L ° A + L 0A C 15 exp(- 7 i 5 L) + C x exp(- 7l L a2 ) 
<C 13 exp(-7 13 L a2 ) 
for appropriate C13 < co and 713 > 0. The lemma follows. □ 

Invading the spatial box Ji{z). By Lemma 3.7, our selected path At 
reaches the target set J$ by time T — \J~L with probability close to 1. The 
last step, similar to the one described in [11], is to show that such a path 
gives us a good opportunity to bring a 1 up to (x,T). 

To define At between time Tjy and time r = T — y/~L, we pretend that 
our selected path follows the path determined by the algorithm of the first 
ancestor starting at (A T .i , T — Tjy). The next step is to prove that, between 

times Tjy and r, the selected path does not escape from J$. We will then let 
the dual process spread out for the remaining yL time units and sprinkle 
the target region J4 where the l's live. 

Lemma 3.8. Assume that a> f3> \ c . Let T = L L5 andr = T - \[~L. It 
then follows that 

P(At $ J 6 for some T^<t<r)< C 16 exp(-7i 6 L ' 2 ) 

for appropriate C\q < 00 and 716 > 0. 

Proof. This is, with some minor modifications, the same as the proof 
of Lemma 3.3. □ 

To complete the proof of Proposition 3.1, we still need to prove that 
(A T ,\^L) will be occupied by a 1, provided J4 is s-good at time 0. Since 
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the selected path has been constructed to block 3's, this 1 will determine 
the type of (x,T) unless another 1 succeeds earlier. This will prove that the 
specialists present in J4 invade the neighboring box J±(z) with probability 
close to 1. 

Lemma 3.9. Assume that a > (3 > A c and that J4 is s-good at time 0. If 
A T E Jq, there are C17 < 00 and 717 > such that 

P((At,Vl) is not occupied by a 1) < Ci7exp(— 717L ' 1 ). 

Proof. We run the dual process starting at (A T ,yL) for yL — L - 2 
units of time. It follows from the Shape Theorem (see [17], page 128) and 
the ergodicity of the upper invariant measure of the contact process (see 
Proposition 2.16, page 143, in [16]) that, except for a probability smaller than 
Ci8 exp(— 718-L ' 1 ), we can select L ' 1 sites at time L ' 2 which are contained 
in the dual process starting at (A T , y/~L ) such that (i) all these sites are 
contained in J5, (ii) they are at least L 0,3 units apart from each other, 
and (iii) none of the duals starting at these sites uses the same part of the 
graphical representation for the remaining L 0,2 units of time. Since each of 
the duals has a positive probability of surviving, and the square J4 is s-good 
at time 0, each of these L 01 sites has probability n > of being occupied by 
a 1. In conclusion, 

P((A T ,VZ) is not occupied by a 1) < Gi 8 exp(-7i 8 L ai ) + (1 - rj) 1 "' 1 . 

This completes the proof of the lemma. □ 

The proof of Proposition 3.1 is now straightforward. First of all, by com- 
bining Lemmas 3.7, 3.8 and 3.9, it follows that if J4 is s-good at time 0, 
then, for any x £ Ji(z), 

P(£t(x) = 3) < C13 exp(-7 13 L a2 ) + C 16 exp(- 7l6 L a2 ) + C 17 exp(-7 17 L - 1 ) 

for some Cig < 00. In particular, since there are (L/2 + l) 2 sites in J&(z), it 
follows that 

P(£r(x) = 3 for some x € J 4 (z)) < CigiT 1 < e/2 

for sufficiently large L. In other respects, the process dominates a one-color 
contact process with parameter j3 > A c , so the probability that there ex- 
ists an L ' 1 x L ' 1 square D(w) C J^(z) that is empty can be bounded by 
C 2 o^ L8 exp(-72 L ' 2 ) < e/2. Proposition 3.1 then follows. 

Proof of Proposition 3.2. The proof of Proposition 3.2 is almost the 
same. The construction of the selected path, however, is slightly different, 
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since there are no s-arrows in the squares B(0, 1) and 5(1,0). We will use 
the same notation to describe the construction of the new selected path. 

Let ei = (1,0) and x £ Ji{e\). As previously, A t will be the first suc- 
cessful path in a sequence of test paths. The first difference is that, since 
B(e\) C Ti.2, and the l's cannot give birth through an arrow that points into 
Ti.2, we can now apply the repositioning algorithm at each renewal point, 
regardless of whether it is associated with an s-arrow or not. On the other 
hand, the s-arrows are now prohibited, so we remove them from the graphi- 
cal representation. We run the dual process, starting at (x,T) for T units of 
time, and break up the path of the distinguished particle at renewal points. 
Since > A c , the event that the dual process lives forever still has a positive 
probability, even after removal of the s-arrows. 

To define the 1st test path we first set cr^o = 0. The process Q starts 
at £q = x and follows the path of the first ancestor of (x,T) until the first 
time, denoted by c^i, that it jumps to a renewal point. As previously, the 
position of the test path at time <7i 5 i is determined by a repositioning al- 
gorithm. The process then follows the path described by the algorithm of 
the first ancestor until the next time, denoted by o"i 5 2, that it jumps to a 
new renewal point, and so on. Let B ai k denote the second ancestor, if it 
exists, of the previous renewal at time cr l fe . Since (3 > A c , there is a positive 
probability that (B ai k ,T — a\^) lives forever. To describe the new reposi- 
tioning algorithm, we define as above an embedded jump process 5/ which 
stays put except at times di^. Let e 2 = (0, 1) and let A' denote the straight 
line going through Lz and 2Le2- The jump process evolves according to 
the following rules. If site B ai k does not exist, or if the dual starting at 

(B ai k ,T — 0"i,jt) does not live forever, then we set = . Otherwise, 
we have the following three alternatives: 

1. If dist(^ k ,A')> 5 and dist(B ai k , A') < dist(f ^ fc , A'), then we set fc = 

2. If dist(^ lfe , A') < 5 and dist(# CTl k , 2Le 2 ) < dist(|^ fc , 2Le 2 ), then we set 

S °i,k = B <n,k- 

3. Otherwise, we set Si , = f i, , . 

In other words, we now move the jump process toward the target set Js( e 2) 
in a direction perpendicular to the straight line A. 

As previously, the 2nd test path follows the path of the first ancestor of 
(x, T) until the 1st test path reaches Jg(e2), at which time it moves according 
to the repositioning algorithm, and so on. We now call a successful path 
if a particle in Jsi^) following this path by going forward in time does not 
cross any arrow that points to 7i\ on its way up to (x,T). In particular, if 
this particle is of type 3, it will determine the type of (x,T) unless another 
3 succeeds earlier. 
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It is not difficult to see that, with this new construction, all our estimates 
still hold. In particular, if the square Ja{^2) is g-good at time 0, then 

P(£ T {x) = 1 for some x £ J 4 (z)) < C 2 \L~ X 

for a suitable C21 < 00. Proposition 3.2 and Theorem 2 then follow. □ 
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